Effect of phase fluctuation and dephasing on the dynamics of entanglement generation 

in a correlated emission laser 
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A detailed study of the eflects of phase fluctuation and dephasing on the dynamics of the entangle- 
ment generated from a coherently pumped correlated emission laser is presented. It is found that the 
time evolution of the entanglement is significantly reliant on the phase fluctuation and dephasing, 
particularly, at early stages of the lasing process. In the absence of external driving radiation, the 
degree of entanglement and intensity turns out to attain a maximum value just before starting to 
exhibit oscillation which dies at longer time scale. However, in case the driving mechanism is on, the 
oscillatory nature disappears due to the additional induced coherent superposition and the degree 
of entanglement would be larger at steady state. Moreover, the degree of entanglement as predicted 
by the logarithmic negativity and the Duan-Giedke-Cirac-ZoUer criteria exhibits a similar nature 
when there is no driving radiation, although such a trend is eroded with increasing strength of the 
pumping radiation at longer time scale. The other important aspect of the phase fluctuation and 
dephasing is the possibility of relaxing the time at which the maximum entanglement is detected. 

PACS numbers: 42.50.Ar, 42.50.Gy, 42.50.Lc, 03.65. Ud 
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I. INTRODUCTION 

In recent years, a nondegenerate three-level laser has 
attracted a great deal of interest in relation to its poten- 
tial as a reliable source of robust entangled light [H-Q. 
The nonclassical features of the radiation are predomi- 
nantely attributed to the atomic coherence induced by 
initially preparing the atoms in a coherent superposition 
of the energy levels between which a direct spontaneous 
transition is electric dipole forbidden or by pumping 
the same with an external radiation [3|-l6|- However, due 
to the fragility of the coherent superposition, the perti- 
nent continuous variable entanglement is believed to be 
significantly ruined by various phenomena closely asso- 
ciated with the lasing process; noteworthy among many 
are dephasing and phase fluctuations. In practical situa- 
tions, it is not possible to prepare the coherent superposi- 
tion with an arbitrary perfection due to instability of the 
laser, faulty setups, and the inefficiency of the devices em- 
ployed [7Hll|. Even, in case one manages to prepare the 
coherent superposition somehow, it will subseque ntly de- 
cay due to the interaction with the surrounding [ij [ll] ■ 
Therefore, a study of the effects of the phase fluctuations 
and dephasing on the evolution of the entanglement gen- 
eration appears to be of a paramount importance. 

Over the years, the issue of phase fluctuations has at- 
tracted a considerable attention [ijj. In these studies 
the phase and photon number operators are presumed to 
fulfill a canonical commutation relation. In a more re- 
cent times, assuming the phase as randomly fluctuating 
c-number parameter, it is shown that the phase fluctu- 
ations substantially decrease the degree of entanglement 
0) EBl • Investigation of a similar nature where the co- 
herence is assumed to be partially prepared is also con- 
sidered (sl-fiot. In the present contribution, in order to 
investigate the effects of the phase fluctuations and de- 



phasing on the time evolution of the generated entangle- 
ment in depth, coherently pumped correlated emission 
laser turns out to be a convenient tool (detailed descrip- 
tion of the physical model can be found, for instance, in 
piij). In this model, the three- level atoms are taken to be 
prepared initially in a partial 50:50 probability to be in 
the upper and lower energy levels. The atoms prepared in 
this manner are injected into the cavity at a constant rate 
and externally driven with a resonant coherent radiation 
while traversing the cavity. It is not difficult to envision 
that the pumping radiation could have also led to phase 
fluctuations; whose effect was thoroughly addressed else- 
where [3,[la[l3l- Moreover, taking the recent advance in 
locking a phase in such a system into consideration, the 
rate at which the coherence superposition decays is set 
to be less than the spontaneous atomic decay rate. 

In the midest of the efforts geared towards understand- 
ing how to create, manipulate, explain and relate entan- 
glement with various nonclassical correlations, it becomes 
evident that the degree of detectable entanglement when 
different criteria are used is found to be different. The 
emerging disparity can be fundamentally linked to the 
difference in the physical contexts and assumptions in- 
volved in deriving the correlations (T8| . In view of this, 
over the years, quite a large volume of work aimed at 
quantif ying the degree of available entanglement has ap- 
peared [l^l- Essentially, the proposed criteria are in the 
form of some sort of inequality, in most instances provide 
only sufficient conditions and on top of that practical re- 
alization is mostly intractable due to the involved im- 
pairment in the corresponding measurement strategies. 
Therefore, it would be imperative studying the nature 
of entanglement applying different criteria with the in- 
tention of looking for the overlap so that the likelihood 
of the prediction and realizeability becomes more reli- 
able [2^. With this conviction, the time evolution of 
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the entanglement is studied applyin g th e criteria usu- 
ally dubbed as logarithmic negativity j2l[, Duan-Giedke- 
Cirac-ZoUer (DGCZ) [H and Hillery-Zubairy (HZ) [13]. 
The obtained results are compared whenever possible. 

In the study of the time evolution of a similar system 
when many parameters are involved, the usual approach 
is numerically calculating the required correlations mak- 
ing use of either the characteristic function _7] or the 
rate equation [13, [ill- However, in this article, the dy- 
namics of the cavity radiation is analyzed following a 
straightforward analytic procedure recently outlined else- 
where The resulting analytical solutions were 
successfully employed earlier in the study of the time 
evolution of the two- mode squeezing and intensity @. 
Particularly, the time evolution of the two-mode squeez- 
ing shows that the phase fluctuation inherent with the 
imperfect preparation can enhance nonclassical features 
when the atoms are driven with a strong external radia- 
tion by creating spontaneous emission transition root via 
weakening the otherwise induced strong coupling of the 
upper and lower energy levels by pumping. However, in 
this contribution, it is shown that the degree of entangle- 
ment exhibits damping oscillatory nature when there is 
no driving radiation where in between there is a narrow 
window of time in which the degree of entanglement is 
found to increase with the deviation of the phase fluctua- 
tion and rate of dephasing. The fluctuation of the degree 
of entanglement is solely attributed to the choice of the 
values of these parameters which are not completely in- 
dependent. 



II. EQUATIONS OF EVOLUTION 

The interaction of a pumped nondegenerate three-level 
cascade atom with a resonant two-mode cavity radiation 
can be described in the rotating- wave approximation and 
the interaction picture by the Hamiltonian of the form 

H = ig[a\a) {b\ - \b) {a\a) + b\b) (c| - |c) {h\P] 

+ ^^[\c){a\-\a){c\], (1) 

where f2 is a real-positive constant proportional to the 
amplitude of the driving radiation and <? is a coupling 
constant chosen to be the same for both transitions, a 
and b are the annihilation operators that represent the 
two cavity modes. In the cascade configuration, the 
transition from upper energy level \a) to the interme- 
diate level \h) and from level \b) to the lower energy level 
|c) are taken to be resonant with the cavity radiation, 
whereas the transition \a) O |c) is electric dipole for- 
bidden. Moreover, as a result of the emerging various 
quantum effects, it is assumed that the atoms can only 
be initially prepared in a partial maximum coherent su- 
perposition of the upper and lower energy levels (a more 
general case has been presented in [1]). 

On account of this, the initial state of the three-level 



atom can be written as 

|*a(0)) ^^[\a) + e^^\c)l (2) 

where tp is an arbitrary phase randomly distributed about 
a fixed mean value (/?o- Based on this assumption, the 
phase can be defined as ip — i^q + dtp m which 5ip is taken 
to be small random fluctuations around tpo that can be 
adjusted at will by proper choice of the phase of the cav- 
ity radiation [lO] . In line with Eq. ([2]) , the corresponding 
initial reduced atomic density operator appears to be 

Pa(0) = \[\a){a\ + e-'^^\a){c\ + e^\c){a\ + \c){c\] , (3) 

where e^"^/2 stand for the initial atomic coherence. 

Addressing the contribution of every phase change does 
not seem practically realistic. Consequently, assu ming 
the phase fluctuations as a Gaussian random process |10l . 
[lll,[ll and using the deviation of the phase fluctuation 
instead of the actual phase led to more tractable situation 
^9]. With similar conviction, employing the property of 
the Gaussian variables [l^l , one can express 

(exp±i(5(p) = exp -((5(/5^/2). (4) 

Moreover, for a Gaussian random process {54>) is zero 
where ((5^0/2) — 9 represents the deviation which is gen- 
erally designated as phase fluctuation. 

Furthermore, following the detailed proceduce of ob- 
taining the master equation and equations of evolution 
presented elsewhere [8, 9], it is possible to verify that 

a{t) = C+(t)«(0) + D+it)P*iO) + E+it) + F+{t), (5) 

p{t)^C-it)(3{0)+D-it)a*iO) + E-{t)+F-it), (6) 
where 

C±(i) = i [(1 ±p)e-^-* + (1 Tp)e-^-^*] , (7) 

D±{t) ^ ^[e-^-' - e-'^-% (8) 

-f (l-p)e-'^+(*-*')]/,(t')di', (9) 

+ (l+p)e-^+(*-*')]/b(t')di', (10) 
^+^^^ " ^ Jj'~'^^^'~''^ ~ e-''- (*-*')]/,* (t')di', (11) 
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F-it) 



-e-f'-^'-'"^]f*{t')dt\ (12) 



with 



±^r-(i + KV + + Xf - [(2 - C'Oe- 
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(13) 



VC'2(1 + CO' + 4(C'2 + X)^ - [(2 - C'Oe-'] 



(14) 



9± 



-C'(l + C'C)±[(2-C'C)e-^] 



VC''(1 + CO' + 4(0 + Xy - [(2 - coe-'v ' 

(15) 



in which 



2ra5' 



i? = (4 + c')(i + C'C), 



(16) 



(17) 



( — ft/'^, C' = ^^/r, and X — 7/r. It is worth noting that 
K is the cavity damping constant, F is the spontaneous 
atomic damping rate, 7 is the rate at which the coher- 
ent superposition decays and is the rate at which the 
atoms are injected into the cavity. 

In addition, the corresponding stochastic noise forces 
are shown to have correlations of the form 

{fait')f:it)) = ^ [20 + 2x - (2C' + C)e-^] Sit - t'), 

(18) 



(A(t')A*W>=o, 



(19) 



A 



{hit') fait)) - — [C'(l + CO + (2 - C'Oe-"] 6it - t'), 

(20) 



{f:{t')Ut)) = ifaitVait)) = {Mt')Mt)) = 0. (21) 



III. CONTINUOUS VARIABLE 
ENTANGLEMENT 

In recent years, on the basis of different conditions and 
assumptions, several inseparability criteria for continuous 
variable product states have been proposed (2^, [2^ [2^ . 



One of these is related to a logarithmic neg ativity (E'at) 
[Hill [13 which can be defined as |27H29| 



En = max[0,-log2 Vs], 



(22) 



where Vs is the smallest eigenvalue of the symplectic ma- 
trix [23, [2^. The inseparability condition associated with 
this approach is based on the well established fact that 
for the composite state to be separable to its constitute, 
the product density operator should have a positive par- 
tial transpose. 

In light of this, upon solving the eigenvalue equations 
for the symplectic spectrum of the covariance matrix of 
the partially transposed density operator, the smallest 
eigenvalue is found to have a form 



^-^e- 4detS 



1/2 



with 



^ = detcTA + detfTB — 2detcrAB , 



(23) 



(24) 



where a a and as are the covariance matrices describing 
each modes separately while ctab represents the inter- 
modal correlations |2l| . 

In connection with the bimodal states involved, the 
2X2 block form of the covariance matrix S can be ex- 
pressed as [S^l 



^AB '^B 



where 



-{X,Xj+X,X,) -{X,){X,). 



(25) 



(26) 



The corresponding quadrature operators are defined as 
Xi = a + at, X2 = i(at - a), X3 ^ b + and X4 = 
iib^ — b). With this introduction, the covariance matrix 
can be expanded as 



where 



m c 
m — c 
c n 
-c n 



^l + 2{a\t)ait)) + {a\t)) + {a^\t)) 
-[(a(t))' + (at(t))2+2(a(i))(at(0)], 

, = l + 2(bHt)bit)) + {b'it)) + {b^\t)) 
-[{bit)f + {bHt))^ + 2{bit)){b\t))], 



(27) 



(28) 



(29) 



: {ait)bit)) + (at(t)&t(t)) + (a(t)6t(t)) + (at(t)S(i)) 

[{ait))m) + {aHt)){bHt)) + m){bHt)) 
{&Ht)){bim- (30) 
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These equations can be expressed in terms of c-number 
variables associated with the normal ordering as 



m = 1 + 2{a*{t)a{t)), 



n-l + 2(r(i)/3W), 



(31) 



(32) 



+ A 
-A 



Lql+Mq_{l+p) 



8Bn+ 
Lql - Mq^{l - p) 

Lq^_ + Mq^p 



[1 - e-^A'+t] 
[1 - e-^A"-*] 



[1-e 



-(p++Ai-)f 



, (40) 



c=(a(0/3(0) + («*(i)r(i)). 



(33) 



In obtaining these expressions, upon neglecting the inter- 
action among the atoms in the cavity and assuming the 
cavity to be initially in a vacuum state, it is taken that 
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y{t)) = m) = {c^\t)) = {P\t)) = {a{t)P*{t)) = 0. 



Furthermore, making use of Eqs. ([25]), (l26l). ([27 | . dST]) . 
32|) and (|33|) . one can readily see that 



detcTA = 1 + ^{a*{t)a{t))[{a*{t)a{t)) + 1], (34) 



detae = 1 + {t)(i{t))W {t)m) + 1], (35) 



detS = 16 



detfjAB = -4(a(t)/3(t))2, 

1 {a*{t)a{t)) + {(3*(t)(i{t)) 
4 2 



(36) 



{a*{t)a{t)){p*{mt)) - Hmt)f 



(37) 



It is a well established fact that the two-mode Gaus- 
sian product state would be entangled if En = — Log2Vs 
which implies that Log2Vs should be negative. That 
means 



Vs < 1 



(38) 



represents the condition for detecting entanglement. In 
order to quantify entanglement using this approach di- 
rectly, the various correlations in Eq. (1371) should be 
determined applying Eqs. ^ and ([6]). To this end, as- 
suming the cavity to be initially in a two-mode vacuum 
state and the noise force at time t does not statistically 
related to the cavity mode variables at earlier times, one 
can readily verify that Q 



{a*{t)a{t))=A 
+ A 
+ A 



L{l-pf+Mq+{l-p) 



8Bfi+ 

Ljl+py^-Mq+jl+p) 

L{l-p^) + Mq+p 
2B(/i++A*-) 



fi- 



ll- 



[1 - e-(A'++P-)t]^ 
(39) 



t{t)m) = A 



A 



2£g_(l-p) + i\/(l-p^+g_g+) 
16B/i+ 

2Lq_{l+p)-M{l-p^ + q,q+Y 



X [1 - e-^A'+t] 



16B/i_ 

X [1 - e-2M-t] 

'2Lq^p + M{l+p^ -q-q+) 



A 



where 



4S(/x+ + /i_) 

L = 2C'2 + 2x-(2C' + C)e"' 
M = C'(l + C'C) + (2-C'C)e- 



(41) 

(42) 
(43) 



In the following, the time evolution of the degree of de- 
tectable entanglement is investigated employing the cri- 
terion of logarithmic negativity (|38l) . To achieve the in- 
tended goal, Vs is plotted against time (a scaled unitless 
parameter) applying Eqs. (Hg, dM]), ([Ml), dM]), 
([37| . ([39]). ffl . dHJ, ^ and (gH). In order to see the 
nature of entanglement in depth, various special cases are 
considered. 



A. When there is no external driving radiation 

In recent years, the investigation of the quantum fea- 
tures of the radiation generated by the nondegenerate 
three-level cascade laser when 7 = F and 9 = Q has at- 
tracted interest. These studies are usually confined to 



the steady state case 
lected parameters [1, 



[31 1 or numerical solution for se- 
However, in this article, upon 




using the rate of dephasing and deviation of the phase 
fluctuation as a probe, the evolution of the entanglement 
is studied when f2 = 0. 

As clearly shown in Figs. [U [21 [3] and lU the degree 
of entanglement exhibits damping oscillatory nature for 
certain parameters under consideration. It is worth not- 
ing that, since the injection of the atoms is launched at 
t = 0, entanglement would not be expected at the begin- 
ning when the cavity is taken to be initially at vacuum 
state. Moreover, it is not difficult to see that the degree 
of entanglement increases rapidly and then starts to fluc- 
tuate with time. The increment of the degree of entangle- 
ment at early stages of the lasing process can be directly 
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FIG. 1: Plots of the time evolution of the smallest eigenvalue 
(Vs) of the cavity radiation for k = 0.5, 6 = 0, A = 10, Q, — 
and different values of 7/r. 



FIG. 3: Plots of the time evolution of the smallest eigenvalue 
(Vs) of the cavity radiation for k — 0.5, 7 = 0.75r, A = 10, 
f2 = and different values of Q. 
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FIG. 2: Plots of the time evolution of the smallest eigenvalue 
(Vs) of the cavity radiation for k = 0.5, Q — 0.25, A — 10, 
f2 = and different values of 7/r. 



linked to the increment of the atoms that participate in 
the correlated emission. Similar increment behavior at 
early stages of lasing operation was observed for various 
schemes of three-level laser under different conditions us- 
ing various entanglement criteria [1, 0, [IS Ell • 

Nevertheless, in order to trace the origin of the fluc- 
tuation, it appears natural referring back to Eq. (1131) . 
In Eq. (|13p . if 17 = 0, the term under square root re- 
duces to \/4:{x^ — e^2^). This indicates that, in Eqs. 
(p9|) . (j40|) and (|4T1) . the exponential term that contains 
time would be complex valued which leads to oscillation 
if 7/r < . In case of = and 7 = F, generally the 
oscillation in the entanglement evolution is not observed 
as large volume of earlier works asserted [1, 0, lol-flll [l5l | . 
One can hence deduce that the main source of this os- 



FIG. 4: Plots of the time evolution of the smallest eigenvalue 
(Vs) of the cavity radiation for k — 0.5, 9 — 0.25, 7 — 0.75r, 
S7 = and different values of A. 

dilatory nature is the competition between the rate at 
which the coherent superposition decays and the degree 
at which the same coherence exhibits fluctuations ini- 
tially. Quite obviously, this competition can be directly 
associated with the mean photon number of the gener- 
ated radiation. Although it is not presented in this work, 
the pertinent mean photon number also exhibits a sim- 
ilar kind of oscillatory nature. Therefore, comparison 
between the two reveals that the degree of entanglement 
and the mean photon number are larger for the same pa- 
rameters. It is hence anticipated that the thermal fluc- 
tuations arising due to high intensity might be responsi- 
ble for the downturn of the quantum features (increment 
of 14) 0. It is worth noting at this juncture that in 
the definition of the covariance matrix, terms related to 
normal ordering are included that are basically associ- 
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ated with the contribution of vacuum fluctuations which 
might have played a significant role. 

It is not hard to notice from Fig. [T] that the small- 
est eigenvalue corresponding to the covariance matrix is 
smaller for smaller values of 7/r at early stages of the 
lasing process. This indicates that the degree of entan- 
glement by large decreases with the increasing rate at 
which the coherence superposition decays provided that 
the rate at which the atoms spontaneously decay is taken 
to be constant. However, it would be appropriate noting 
that for small pockets of time, the degree of entanglement 
can increase with the rate at which the coherence decays. 
For relatively longer time span, the degree of entangle- 
ment slowly oscillates very close to Vg = 0.5. This can 
be fairly related to the earlier report for the same system 
when 7 = r, ?7 = 0, = and A = 10, where Vg is found 
to be 0.5 at steady state [l^l- On the other hand, critical 
scrutiny of Figs [T] and [5] reveals that the time at which 
the degree of entanglement becomes maximum increases 
with the rate of dephasing. Though increasing the rate of 
dephasing significantly damages the generated entangle- 
ment, it can be utilized in prolonging the time at which 
the maximum entanglement can be generated. This can 
be taken as a positive aspect in practical realization of 
entanglement in such a system where the corresponding 
actual time is very small. 

In plotting Fig. [51 two changes are made; the val- 
ues of 7/r are reduced while the value of 9 is increased 
so that the nature of the evolution of entanglement be 
clearly evident. If one expects the nature of entangle- 
ment to follow the same trend for various values of 7/r, 
the degree of entanglement indicated in Fig. [2] should 
have been much larger. However, when the results pre- 
sented in Figs. [2 and [5] are compared, the entanglement 
turns out to be better in the former. This entails that 
the phase fluctuation significantly ruins the obtainable 
degree of entanglement at earlier times. For relatively 
longer time span, the oscillation is limited to very close 
to Vs = 0.5 as in the case when 9 — 0. This suggests that 
the effect of phase fluctuation can be minimal at longer 
time scale when compared to earlier times. In order to 
attest to this claim, Vs is plotted against time for differ- 
ent 9. The result depicted in Fig [3] indicates that, indeed, 
the degree of entanglement decreases with the deviation 
of the phase fluctuation in the earlier times of the lasing 
operation. In relation to the involved oscillation, there 
is also a brief window of time in which the degree of en- 
tanglement can increase with the phase fluctuation. Just 
as we have observed for the rate of dephasing, the time 
at which the maximum entanglement can be detected in- 
creases with the phase fluctuation. 

Aiming at the prospect of increasing the degree of en- 
tanglement, Vs is plotted against time for 9 = 0.25, 
7 — 0.75r and different values of the linear gain coef- 
ficient. As clearly shown in Fig. |4l the degree of entan- 
glement increases with the linear gain coefficient in the 
early stages of the lasing process. Quite remarkably, in 
light of the oscillation, there is a window of time where 



the entanglement can decrease with the linear gain coeffi- 
cient. Once again, at larger time scale the dependence of 
the entanglement on the linear gain coefficient tends to 
be much lesser when compared to the earlier times. Fit- 
ting things together may indicate that the nonclassical 
features and intensity of the cavity radiation take max- 
imum values before they start damping oscillation that 
dies at longer time scale where the steady state properties 
begins to take effect. 



B. In stronger driving regime 

The entanglement generated from coherently pumped 
correlated emission laser when 77 = 0, 7 = F and 9 — 
was studied earlier using the criteria of DGCZ. It is gen- 
erally noted that Au^ + Aw^ is very close to one at steady 
state in a strong driving limit (fi >> F) [l^. In the fol- 
lowing, upon varying the values of the rate of dephasing 
(7/r) and the deviation of the phase fluctuation {9), the 
time evolution of the entanglement that can be detected 
by the logarithmic negativity is studied. 




y/r=0.80 

y/r=0.85 

y/r=0.90 



1 2 3 4 5 

time 

FIG. 5: Plots of the time evolution of the smallest eigenvalue 
(Vs) of the cavity radiation for k = 0.5, 9 = 0, A = 10, 
fl = lOr and different values of 7/r. 

In line with earlier discussion, the term under square 
root in Eq. ([T^ would not be negative for large il. As a 
result, as clearly evinced in Figs. [5l|6l[7]and[8l the degree 
of entanglement does not exhibit oscillatory nature in the 
strong driving limit. It is, rather, observed that the de- 
gree of entanglement increases with time for all param- 
eters under consideration. This outcome suggests that 
a better degree of entanglement would be expected at 
steady state which is believed to be an important aspect 
in practical utilization of the system. It can readily be 
seen from Fig. [5] that when the atoms are pumped exter- 
nally with strong coherent radiation, the entanglement 
appears to decrease with the rate of dephasing. More- 
over, one can observe from Fig. [5] that the nature of the 
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FIG. 6: Plots of the time evolution of the smallest eigenvalue 
{Vs) of the cavity radiation for k = 0.5, 9 = 0.25, A = 10, 
fl = lOr and different values of 7/r. 
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FIG. 7: Plots of the time evolution of the smallest eigenvalue 
(Vs) of the cavity radiation for « = 0.5, 7 = F, y4 = 10, 
Q, = lOF and different values of 9. 



degree of entanglement does not display much change 
from what is imparted in Fig. [5] except decrement in the 
effect of the rate of phase fluctuation slightly. In order 
to show the dependence of the evolution of the entan- 
glement on the phase fluctuation in depth, Vs is plotted 
against time for different values of 0. Relatively spaced 
values of 6 are chosen so that the modification in the 
evolution of the entanglement be evident from the plots. 
As can easily be noticeable from Fig. [3 the degree of 
entanglement increases with the phase fluctuation. This 
is in a complete agreement with a recent observation that 
in case when there is a strong external driving the phase 
fluctuation enhances nonclassical features via creation of 
indirect spontaneous emission root 

It is not difficult to see from Figs. [SI [Hand [7] that Vs is 
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FIG. 8: Plots of the time evolution of the smallest eigenvalue 
(Vs) of the cavity radiation for k = 0.5, 7 = 0.75r, 9 = 0.25, 
fl = lOF and different values of A. 



very close to 0.4, which is by far larger than the case when 
= 0. In the quest of improving the degree of entangle- 
ment, Vs is plotted against time for 7 = 0.75F, 9 = 0.25 
and different values oi A. As unambiguously presented 
in Fig. m the degree of entanglement is found to increase 
with the linear gain coefficient. This implies that in the 
process of practical utilization, one can manipulate the 
rate at which the atoms are injected into the cavity to get 
an optimum entanglement. It is worth noting that the 
values of A are set in this manner so that one can able to 
make comparison with the case when i7 = rather than 
aiming at finding the highest degree of entanglement. In 
light of this, upon comparing the results shown in Figs, 
m and m one can readily see that the degree of entangle- 
ment is larger when = at earlier times, whereas the 
opposite holds for longer time span where the oscillation 
naturally subsides. Although there is no earlier analysis 
that supports this outcome, the result obtained here en- 
tails that the driving mechanism is required if one wishes 
to generate strongly entangled light at steady state for 
the system under consideration. 



IV. COMPARISON OF THE DEGREE OF 
ENTANGLEMENT 

It is a well known fact that the criteria used to de- 
tect the degree of entanglement of the continuous vari- 
able composite system are of different kind and also re- 
lated to different physical contexts. It so happens that 
generally they do not lead to similar conclusion since pre- 
dominantly they are emanated from different principles 
and assumptions. In order to get a clear understand- 
ing of the nature of entanglement, it would be necessary 
studying the relation among the predictions of various 
criteria. With this understanding, in the following, the 
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evolution of the entanglement that was obtained apply- 
ing the criterion of logarithmic negativity is compared 
with the DGCZ and HZ criteria. 



A. Duan-Giedke-Cirac-Zoller criteria 

In this section, the nature of entanglement as pre- 
dictable by the criterion following from the variance of 
the Einstein-Podolsky-Rosen (EPR) type quadrature op- 
erators [22] is compared with logarithmic negativity that 
is employed in the preceding section of this work. Pre- 
vious study shows that the nature of the predicted de- 
gree of entanglement by the two criteria for a similar 
system when fl — 0, 9 = and 7 = F has the same form 
except near rj = where significant disparity was ob- 
served at steady state [2^ . The main aim here is shifted 
to the study of the time evolution of the entanglement 
upon using 7/r and 9 as a, probe. It is worth noting 
that the detailed analysis of the degree of entanglement 
that can be quantified by the criterion following from 
DGCZ has been presented in connection with the degree 
of two-mode squeezing elsewhere @ . In earlier work, the 
oscillatory nature of the entanglement was not observed 
since the amplitude of the driving radiation was set to be 
il > O.Sr. With the aid of the calculation presented ear- 
lier, the variance of the quadrature operators (Au^-t-Aw^) 
and the smallest eigenvalue of the symplectic matrix (Vs) 
are plotted against time for various parameters. In the 
plot {Au^ + Av'^) is multiplied by half so that comparison 
can be as straightforward as possible. 




FIG. 9: Plots of the time evolution of the smallest eigenvalue 
(Vs) and variance of the EPR operators (Au^ + An^) of the 
cavity radiation for k — 0.5, 6 = 0.25, A — 25, Q — and 
7 = 0.75r. 

It is unequivocally shown in Fig. [9] (please note that 
A = 25 is used so that the dependence of the two can be 
distinctly observed) that the degree of entanglement that 
can be predicted by the two criteria exhibits very similar 
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FIG. 10: Plots of the time evolution of the smallest eigenvalue 
(K) and variance of the EPR operators (Au^ -f An^) of the 
cavity radiation for k = 0.5, Q = 0.25, A = 25, Vl = lOF and 
7 = 0.75r. 



nature when f2 = 0. However, their nature turns out to 
be completely different when 57 7^ in the longer time 
span as evinced in Fig. [TU] where the time parameter is 
limited to small interval so that the relation between the 
two can be readily noticeable. It is worth noting that 
Au^ -|- Av^ rapidly increases with time which leads to 
dying of entanglement whereas Vg continuously decreases 
with time. What can be said at this juncture is that the 
two criteria predict similar nature of entanglement evo- 
lution for relatively smaller time scale, but they can also 
predict entirely different thing when there is an external 
driving radiation at longer time span. By the way, such 
a disparity between the two is not uncommon for 7^ = 
which includes the present case [20j . 



B. Hillery-Zubairy criteria 

In earlier communication, on account of the relation 
of the criterion introduced by Hillery and Zubairy 23 1 
with the correlation of the photon number, the gener- 
ated entanglement has been tipped to be quantifiable 
via simultaneous two-photon count measurement [l2| . In 
search of more appropriate and easy means of detecting 
entanglement, the degree of entanglement that can be 
predicted by the logarithmic negativity and HZ criteria 
has been compared earlier for nondegenerate three-level 
laser when = 0, = and 7 = F. It was shown that 
both criteria predict the manifestation of entanglement 
at steady state, even though the degree of nonclassical- 
ity cannot be compared since HZ-criteria has no upper 
limit. In this contribution, upon changing the above pa- 
rameters, the comparative evolution of the entanglement 
is investigated. 

According to HZ criteria, one can say that there is 
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entanglement if 



\{amm>\/{Na{t)){m)), m 

where {Na{t)) = {a1{t)a{t)) and {Nb{t)) = {b''{t)b{t)) 
are the pertinent photon numbers corresponding to the 
involved modes, and {a{t)b{t)) is the intermodal correla- 
tion. Since the operators are already put in the normal 
order, Eq. (|44p can be expressed in terms of c-number 
variables associated with the normal ordering as 

{a{t)m)' > {a*it)a{t)){l3*{mt))- (45) 

For nondegenerate three-level cascade laser, based on 
the assumption that there is no interaction among in- 
jected atoms, the photon number correlation is found to 
have the form [l| 



9{na,nb) = 1 + 



\a*(t)a{t)){P*{t)m)' 



(46) 



Hence upon comparing Eqs. (|45|) and (|46|. one can read- 
ily see that the HZ inseparability criterion can be restated 
as 



(47) 



It is obvious to see from Eq. (|47p that g{na, rib) > 2 for 
all parameters under consideration. It is also clearly in- 
dicated that the evolution of the generated entanglement 
exhibits oscillatory nature. At this juncture, it is essen- 
tial to note that for earlier time span gijia, Ub) goes to in- 
finity in connection with the mean photon number in one 
of the modes rapidly close to zero. In earlier report, sim- 
ilar behavior of g{na,nb) has been explained in relation 
to the correlation involved in the criteria [23] . In light of 
this, it is possible to infer that except near the beginning 
of lasing operation, the two criteria unequivocally assert 
that the system under consideration generates entangled 
light. The overlap between these predictions may be uti- 
lized in detecting entanglement using the simultaneous 
photon count measurement in case the usual homodyne 
detection is found to be practically challenging. 



V. CONCLUSION 

Detailed analysis of the time evolution of the entangle- 
ment generated from a correlated emission laser is pre- 
sented. Aiming at searching for a more reliable outcome, 
the entanglement criteria associated with the positivity 
of the transpose of the density matrix, the variance of 
the EPR type quadrature operators and photon number 
correlation are employed. Based on the assumption that 
various processes affect the preparation of the coherent 
superpositions and subsequent dynamics, the degree by 
which the phase fiuctuates and the rate at which the co- 
herence decays are taken as the probing mechanism. It 
is generally observed that in the absence of an external 



driving radiation (Q. = 0), when the phase can be locked 
(7 < r) and/or due to imperfect preparation when there 
is phase fluctuations (0 7^ 0) , the generated entanglement 
exhibits some sort of damping oscillatory nature. Upon 
varying the rate at which the coherent superposition de- 
cays and the deviation of the phase fluctuation, the de- 
gree of entanglement is found to be altered significantly 
particularly in the early stages of the lasing process. 

However, for strong external driving radiation, the os- 
cillatory nature of the quantum features and statistical 
properties of the radiation disappear. This is mainly 
because, the strong external driving radiation induces 
additional coherence while the atoms traverse the cav- 
ity which subsequently compensates for the effect of de- 
phasing and phase fluctuations. This explanation can be 
vividly evident if one compares the effect of the two when 

= and Q, ~ lOF. Moreover, the oscillatory nature of 
the nonclassical features and photon statistics disappears 
for longer time span, since the steady state effects start 
to dominate. 

Critical scrutiny reveals that physically valid result can 
be obtained by carefuU selection of these parameters. 
This can be related to the fact that when entirely viewed 
from physical perspective, the rate of dephasing and the 
phase fluctuation cannot be set arbitrarily since these 
physical processes are not entirely independent. This 
is mainly because, setting 7/r small principally implies 
that the rate at which the coherent superposition decays 
is small. This, in other words, means that the phase be- 
tween the upper and lower energy levels is locked by some 
degree. This does not arbitrarily go with the assumption 
that the phase fluctuates vigorously (large phase fluctua- 
tions). What possibly can be done in practical utilization 
of the system is adjusting the two parameters in such 
away that one assumption does not preclude the effect of 
the other. 

Comparison among the prediction of the nature of the 
evolution of the entanglement by various criteria more or 
less indicates the same thing. At this juncture, it would 
be appropriate citing the disparity between the predic- 
tion following from the criteria of logarithmic negativity 
and DGCZ for large at longer time scale. By and 
large, this study evinces that the system under consider- 
ation can be a source of strongly entangled light for wide 
selection of parameters. Even though detailed analysis 
at steady state is still lacking, it is possible to infer that 
at longer time scale, externally driving the cavity signifl- 
cantly improves the degree of nonclassical features of the 
generated radiation. 
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